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Abstract 

We review the derivation of the Gaudin model with integrable boundaries. 
Starting from the non-symmetric R-matrix of the inhomogeneous spin-Vi XXZ chain 
and generic solutions of the reflection equation and the dual reflection equation, the 
corresponding Gaudin Hamiltonians with boundary terms are calculated. An alter- 
native derivation based on the so-called classical reflection equation is discussed. 



*E-mail address: nantonio@math.ist.utl.pt 
+ E-mail address: nmanoj@ualg.pt 
iE-mail address: zoltan.nagy@m4x.org 



I INTRODUCTION 



I Introduction 

Gaudin models have applications in many areas of modern physics, from quantum 
optics HI 13 to physics of metallic nano-grains, see |3| and reference therein. A model of 
interacting spins in a chain was first considered by Gaudin HUHJ. In his approach, these 
models were introduced as a quasi-classical limit of the integrable quantum chains. 
Moreover, the Gaudin models were extended to any simple Lie algebra, with arbitrary 
irreducible representation at each site of the chain HQ. 

The rational s£(2) invariant model was studied in the framework of the quantum in- 
verse scattering method |6|. The quantum inverse scattering method (QISM) l7ll8ll9l lT0ll 
as an approach to construct and solve quantum integrable systems has lead to the the- 
ory of quantum groups ||TTl[T2| . In his studies, Sklyanin used the s£(2) invariant clas- 
sical r-matrix ]6|. A generalization of these results to all cases when skew-symmetric 
r-matrix satisfies the classical Yang-Baxter equation lfl3H was relatively straightforward 
fHHl5ll . Therefore, considerable attention has been devoted to Gaudin models corre- 
sponding to the the classical r-matrices of simple Lie algebras Hl6l W7\ ITBI and Lie su- 
peralgebras IIT9ll20ll2Tll22ll23l . The spectrum and corresponding eigenfunctions were 
found using different methods such as coordinate and algebraic Bethe ansatz, sepa- 
rated variables, etc. Correlation functions for Gaudin models were explicitly calculated 
as combinatorial expressions obtained from the Bethe ansatz [5|. In the case of the s£(2) 
Gaudin system, its relation to Knizhnik-Zamolodchikov equation of conformal field 
theory fl24l l25l l26l or the method of Gauss factorization [27J, provided alternative ap- 
proaches to computation of correlation functions. The non-unitary r-matrices and the 
corresponding Gaudin models have been studied recently, see 11281 and the references 
therein. 

A way to introduce non-periodic boundary conditions compatible with the integra- 
bility of the quantum systems solvable by the QISM was developed in [29 1. The bound- 
ary conditions at the left and right sites of the system are expressed in the left and right 
reflection matrices. The compatibility condition between the bulk and the boundary of 
the system takes the form of the so-called reflection equation |30l[3TH . The compatibility 
at the right site of the model is expressed by the dual reflection equation. The matrix 
form of the exchange relations between the entries of the Sklyanin monodromy matrix 
are analogous to the reflection equation. Together with the dual reflection equation they 
yield the commutativity of the open transfer matrix |29l[32l[33ll34l . 

The starting point to obtain the Gaudin model in the framework of the QISM is 
the monodromy matrix of the corresponding inhomogeneous spin chain [6|. Hikami, 
Kulish and Wadati showed that the semi-classical expansion of the transfer matrix of 
the periodic chain, calculated at the special values of the spectral parameter, yields the 
Gaudin Hamiltonians l35l l36l . Hikami generalized this approach to the case of non- 
periodic boundary conditions l37ll . In his work Hikami used the Sklyanin monodromy 
matrix of the open inhomogeneous spin chain. By the semi-classical expansion of the 
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open transfer matrix, at the special values of the spectral parameter, the corresponding 
Gaudin Hamiltonians with boundary terms were obtained |37]. These results were later 
extended to non-diagonal reflection matrices [38, 39] and generalized for other classes 
of Gaudin models, like for example elliptic ones EDI , as well as to other simple Lie 
algebras [41 J. Analogous results were obtained for the open Gaudin models based on 
Lie superalgebras [42 J. 

As it was noticed above, in the periodic case, the Gaudin models are based on clas- 
sical r-matrices which have a unitarity property. Namely, the density of the Gaudin 
Hamiltonians coincides with the classical r-matrices and the condition of their com- 
mutativity is nothing else but the classical Yang-Baxter equation fl5l H7| . In the case 
of non-periodic boundary conditions the derivation of the generating function of the 
corresponding Gaudin Hamiltonians is still an open problem. Sklyanin has developed 
an approach based on the classical reflection equation, involving a unitary classical 
r-matrix and corresponding reflection matrix |[43ll44| . But the derivation of the corre- 
sponding Hamiltonians is not straightforward. A method based on a non-unitary clas- 
sical r-matrix, obtained from a unitary solution of the classical Yang-Baxter equation 
and the corresponding reflection matrix, was proposed |28l[45l[46l . However, the cor- 
responding Hamiltonians do not coincide with the ones obtained by the semi-classical 
expansion. 

This paper is organized as follows. In Section II the R-matrix of the XXZ model 
and its properties are reviewed. The equivalence between the reflection equation and 
the dual reflection equation is shown and the corresponding reflection matrix is given 
explicitly. In Section III an inhomogeneous XXZ spin chain with N sites is studied. It is 
shown how the semi-classical expansion of the transfer matrix, for special values of the 
spectral parameter, yields the Gaudin Hamiltonians, both in the case of periodic and 
non-periodic boundary conditions. These Hamiltonians are calculated in Section IV. In 
particular, the Gaudin Hamiltonians with the boundary terms are given explicitly. An 
approach to study open Gaudin model based on classical r-matrix, classical reflection 
equation and corresponding reflection matrix is discussed in Section V. 



II Reflection equation 



In the framework of the QISM JZHHJ- 13 HDJ the starting point in the study of a quantum 
solvable system is an R-matrix. Here we consider the R-matrix of the spin-% XXZ chain 

BZHMIH 

I \ \ 



R(A,r,) 








sinh(A) — e~ A sinh(?/) 
sinh(A — Y]) sinh(A — r\) 
— e A sinh(//) sinh(A) 



sinh(A — Vj) sinh(A — tj) 
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here A is the so-called spectral parameter and // is the semi-classical parameter. The 
R-matrix satisfies the Yang-Baxter equation in the space C 2 C 2 <g> C 2 

R 12 (A - ii)R 13 (A)R 23 ( H ) = R 23 { F )R 13 (A)R 12 {A - p), (11.2) 

we suppress the dependence on the semi-classical parameter rj and use the standard 
notation of the QISM [ZHSUHJ. [TUJl to denote spaces Vj,j = 1,2,3 on which corresponding 
R-matrices R^, ij = 12, 13,23 act non-trivially. In the present case V\ = V 2 = V 3 = C 2 . 

This form of the R-matrix is related with the symmetric one R l2 (A,fj) = Ri 2 (A,Tj) 
by the similarity transformation 

R n (Kn) -►Adexp(AS?) J R 12 (A,>/), (II.3) 

with S z = diag(l/2, -1/2) (50H5TH52H53H51. The transformed R-matrix still obeys the 
Yang-Baxter equation due to the U(l) symmetry of the initial R-matrix 

[Sf + S| /J Ri 2 (A / J7)]=0. (II.4) 

Sometime ago it was observed |6] that the Gaudin models are related to the classical 
r-matrix and therefore it is essential that the R-matrix has the semi-classical property 

R{KV) = l + yr(A) + 0(ii 2 ), (H.5) 
here r(A) is the corresponding classical r-matrix 



r(A) 
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(II.6) 



which has the unitarity property 

r 21 (_A) = -r 12 (A), (11.7) 
and satisfies the classical Yang-Baxter equation 

[ris (A), r 23 &)] + [r 12 (A - }i),r 13 (A) + r 23 (^)] = 0. (II.8) 

Moreover, for the purpose of deriving the Gaudin Hamiltonians, it is necessary that the 
R-matrix (|II.1|) is normalized so that 



R(0, V )=V, (11.9) 



where V is the permutation matrix in C 2 <8> C 2 . 
The R-matrix <|II.1|) has the unitarity property 



R 12 (A)R 21 (-A) = 1, (11.10) 
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the PT-symmetry 

R[ 2 (\) = R 21 (A), (11.11) 
and the following crossing symmetry property 

^i^V 1 *'^^' (im) 

where £ 2 denotes the transpose in the second space and the matrix J is given by 

' = o )• < II13 > 

We observe that the matrix 

M = jtj =( e e°v) (IU4) 

commutes with the R-matrix 

[M®M,R(A)] = 0. (11.15) 

Notice that the unitarity property (III. 101) and the crossing symmetry ( |II.12| ) imply the 
following useful identity ll50~ll53H 

R; 2 2 (A)M 1 - 1 r1 1 2 (-A + 2^)Mi =p{A-r ] ,r ] )l, (11.16) 

with 



sinh 2 (A) — sinh 2 (//) 
sinh 2 (A) 

A way to introduce non-periodic boundary conditions which are compatible with 
the integrability of the bulk model, was developed in [29]. Boundary conditions on 
the left and right sites of the system are encoded in the left and right reflection ma- 
trices K~ and K + . The compatibility condition between the bulk and the boundary of 
the system takes the form of the so-called reflection equation ll30~ll3Tl . It is written in 
the following form for the left reflection matrix acting on the space C 2 at the first site 
R-(A) G End(C 2 ) 

R 12 (A - F )K- (A)R 21 (A + F )K 2 (p) = K 2 (fi)R 12 (A + y)iq (A)R 21 (A - ji). (11.18) 

In complete generality the compatibility on the right end of the model is encoded 
in the following dual reflection equation |]29ll32ll33ll34H55ll 



A 12 (A - 1 i)K+ t {X)B 12 {\ + H)X 2 +t ( H ) = K+^)C 12 (A + } i)K+ t {X)D l2 {\ - H ). (11.19) 
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where the matrices A, B, C, D are obtained from the R-matrix of the reflection equation 
(111.181 ) in the following way 

Ai 2 (A) = (Ri 2 (A) t12 ) 1 = D 21 (A), (11.20) 
B 12 (A) = ( feUA))^ 2 = C 21 (A). (11.21) 



However, due to the property (III.16D the dual reflection equation (III. 19b can be written 
in the equivalent form 

R 12 (-A + F )X+(A)M 2 R 2 i(-A l^M^K+iv) = 

K+(p)M 1 R 12 (-A -n- 2//)M 1 - 1 X+(A)R 21 (-A + y). (11.22) 
One can then verify that the mapping 

K+(A) = K~{-\-ij) M (11.23) 

is a bijection between solutions of the reflection equation and the dual reflection equa- 
tion. After substitution of (111.231) into the dual reflection equation (|II22|) and using the 
symmetry property (III. 151) one gets the reflection equation ( 1II.18I) with shifted argu- 
ments. 

Although the K-matrix for the XXZ model is well known 11561 l57l l58l , the solution 
corresponding to the R-matrix <|II.1)> differs from that matrix (see 



K W-[ (e 2A_ e - 2 A )c / + e 2A fl ) (H.24) 

which we normalize 

K-(A) = -Lk-(A), (H.25) 

where 



d(A) = f + a cosh(2A) + (Va 2 + Abe) sinh(2A) / (11.26) 

so that 

K-(A)K-(-A) = 1. (11.27) 
Finally, we observe that matrix K~ (A) does not depend on the semi-classical parameter 

^W=0. (11.28) 

Of] 

The right reflection matrix K + (A) is obtained by substituting dll.251) into dll.231) . It is 
important to notice that 

( limX + (A) ) K~(\) = 1. (11.29) 

In general, the normalization conditions (|II.27|) and ( III.29D are not essential in the study 
of the open spin chain. However, together with (III.5b and dll.91) they enable the semi- 
classical expansion of the transfer matrix which yields the open Gaudin model. 
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III Inhomogeneous XXZ spin chain 

In this section we study an inhomogeneous XXZ spin chain with N sites, characterized 
by the local space Vj = C 2 and inhomogeneous parameter &j. For simplicity we start 
by considering periodic boundary conditions. The Hilbert space of the system is 

H = § Vj = (C 2 ) m . 

In the QISM HZHTQUHIll the so-called monodromy matrix 

T(A) = R m {\ - oc N ) ■ ■ • R i(A - «i) (IH.l) 

is used to describe the system. For simplicity we have omitted the dependence on 
the semi-classical parameter t] and the inhomogeneous parameters {a,,/ = 1, ■ ■ ■ , N}. 
Notice that T(A) is a two-by-two matrix in the auxiliary space Vo = C 2 , whose entries 
are operators acting in T-i. Due to the Yang-Baxter equation (|II2|) . it is straightforward 
to check that the monodromy matrix satisfies the RTT-relations |7j[T0l[9l 

R 12 (A - }i)T(\)T( } i) = T( H )T(\)R 12 (\ - ¥ ). (III.2) 

The above equation is written in the tensor product of the auxiliary space Vq®Vq = 
C 2 (8) C 2 , using the standard notation of the QISM and suppressing the dependence on 
the semi-classical parameter tj and the inhomogeneous parameters = 1, . . . , N}. 

The periodic boundary conditions and the RTT-relations (|III.2|) imply that the trans- 
fer matrix 

£(A) = troT(A), (HI.3) 
at different values of the spectral parameter commute, 

[t(A),*oo] = 0/ (in.4) 

here we have omitted the nonessential arguments. 

Following references [35, 36 1 we observe that, due to the normalization of the R- 
matrix (111.91) , the semi-classical expansion of the transfer matrix, for the special values 
of the spectral parameter, is given by 

Z k = f(A = Uk ) = l + f]H k + 0(?/ 2 ), (III.5) 

where H k are the Gaudin Hamiltonians, in the periodic case. The commutativity of the 
Gaudin Hamiltonians, 

[H k ,H l ]=0, (III.6) 

is ensured by the commutativity of the transfer matrix (|III.4|) and the fact the first term 
in the above expansion is the identity matrix, due to (|II.5I) and (III.9D . 
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In order to construct integrable spin chains with non-periodic boundary condition, 
one has to use the Sklyanin formalism [29|. The corresponding monodromy matrix 
T(A) consists of the two matrices T( A) dill.!} and a reflection matrix K~ (A) (111.251) , 

T(A) = T(A)K-(A)T- 1 (-A) 

= Ron (A -«*)••• R i(A - Xi)K£{\)R 10 (\ + «i) • • • R N0 (A + « N ), (IH.7) 

where, for simplicity we have suppressed the dependence on the other parameters. By 
construction, the exchange relations of the monodromy matrix 7~(A) in Vq (g) Vq are 

R 12 (A - p)T(A)R 21 (A + p)7» = T(/y)Ri 2 (A + p)T(A)R 21 (A - y), (III.8) 

using the notation of [29]. The open chain transfer matrix is given by the trace of T(A) 
over the auxiliary space Vo with an extra reflection matrix K + {\) |29|, 

f(A) =tr (K + (A)T(A)). (III.9) 

The reflection matrix K + (A) (|II.23I > is the corresponding solution of the dual reflection 
equation (|II.22|) . The commutativity of the transfer matrix for different values of the 
spectral parameter 

[t(X),t(n))=0, (111.10) 

is guaranteed by the dual reflection equation (|II22|) and the exchange relations flIII.81 ) of 
the monodromy matrix T(A). 

Analogously to the periodic case, the semi-classical expansion of the transfer matrix, 
for special values of the spectral parameter, yields the Gaudin Hamiltonians with the 
boundary terms |38| 

Z k = f (A = «*) = ! + rjH k + C(?/ 2 ), (111.11) 

where are the corresponding Gaudin Hamiltonians. As in the periodic case, the 
commutativity of the Hamiltonians H k is guaranteed by the equation dlll.lOD and the 
normalization conditions (111.51) , dll.9l) and (M.29D . The Gaudin Hamiltonians will be cal- 
culated in the following section. 



IV Trigonometric s£(2) Gaudin Hamiltonians with boundary terms 

In this section we calculate explicitly the Gaudin Hamiltonians. As shown in 1 351136] , in 
the periodic case, it is straightforward to calculate the first two terms in the expansion 

dLSl) 

Z]c\ v =o = tr (RonK -%)'•' Ro/cK - • • • RoiOfc - oil)) |)7=o 

= tr (Ron(«Ic - «n) • • • Pok • ■ ■ Roi(l>C k - Ot-i)) \r, =0 

= tr (P ok ) = h (IV.l) 
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and the Gaudin Hamiltonians in the periodic case, 



7/=0 



H k = ^\ v =o = E tr o ( r onK -«//)••« 9R °' ( ^ ^ ■ ■ ■ V ok ■ ■ ■ R i(% - ai)) 

+ E tr o ^on(«a: - &n) ■ ■ -V ok - ■ ■ 9Rq/ ^^ — — • ■ • Roi(cck - |,/=0 

= Etr [ + gtro (Po* g- l^oj 

= Er H (« fc -«/)- ( IV - 2 ) 

In the case of non-periodic boundary, we consider the left and right reflection ma- 
trices K~ (A) and K + (\) given by ( III.25D and ( III.23D , respectively. In order to obtain the 
expansion (|III.11[) in ( 1IIL9I ) we specify A = oc k and calculate the first term 

Zfc|?7=o = tr (K + (a. k )T{a. k )) \ v=0 = tr (K^(a k )R 0N (a k - a N ) ■ ■ ■ V ok ■ ■ ■ Roi(«fe - «i) x 
xK Q (tt k )R w (tx k + «i) • • • R k0 {2a k ) ■ ■ ■ R m {oi k + oc N )) |^ =0 

= tr (K -(-a fc )7> fcM a *)) = = 1. (IV.3) 

In the last step above we have used the normalization pi.27|) . The Gaudin Hamiltonians 
with boundary terms are related to the second term in the expansion (|III.11|) . 

H k = ^|„=o = tr (^IhWW) + *<> ^WoI^^Ih 
+ X>o (Ko-(-^)^ 9Ro ' ( ^- ffi;) |^ K -K)) 

+ E tr o k(-%)^ (^) 9R ' o( ! t+g,) ir°) • ( IV - 4 ) 

l^k V "7 / 

In the derivation above we have used the fact that the left reflection matrix K~ ( A) does 
not depend on the semi-classical parameter t] (III.28D . Finally the Gaudin Hamiltonians 
can be expressed in a more concise form 

H k = T k (a k ) + r kMk ~ *z) + J2 K k( a khk(^k + «/)*£"(-«*)/ (IV-5) 
where r«(A) is the corresponding classical r-matrix dll.61 ) and 

r,K) = (^M| J7=0 + tro (K -(-a fc )P ^ M (2« fc )) j . (IV.6) 



V CLASSICAL REFLECTION EQUATION AND OPEN GA UDIN MODEL 



Notice that the second term on the righthand side of the (|IV.5I ) coincides with the 
Gaudin Hamiltonians (|IV.2|) and that the first and the third term are the boundary terms 
depending on the reflection matrices K~ and K + . The algebraic Bethe ansatz is used to 
calculate the spectra of these Hamiltonians l38ll4T| . 



V Classical reflection equation and open Gaudin model 



In this section we discuss an approach to study open Gaudin model based on a classical 
r-matrix, classical reflection equation and corresponding reflection matrix. A classical 
reflection equation can be obtained by substituting dIL5D into ( III.18D , taking into account 
(III28I) , and comparing the terms of the first order in rj. However, in the original notation 
of Sklyanin, using the symmetric r-matrix, the formulas are somewhat more compact 
||43l |44| . To this end we apply the similarity transformation dll.31 ) on the classical r- 



matrix dll.611 and obtain the symmetric r-matrix 



r(A) 





I ° 
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1 





cosh(A) 


-1 





sinh(A) 





-1 


cosh (A) 







V o 








o / 



(V.l) 



The unitarity conditions now reads 



r(-A) = -r(A). 



(V.2) 



This r-matrix satisfies the classical Yang-Baxter equation (III.8I) . Together with the corre- 
sponding K-matrix 



K(A) 



d(A) V (e 



e~ A a + e A f (e 



2A 



„2A 



e- 2A )c 



-2A 

e A a + e~ A f 



)b 



(V.3) 



the function d(X) given in dll.26l l, it satisfies the classical reflection equation obtained 
by Sklyanin Il43ll44l 

[r 12 (A - p),Ki(A)K2(fO] + Ki(A)n 2 (A + f/)K 2 (^) - K 2 (^)r 12 (A + ^)Ki(A) = 0. (V.4) 

It is straightforward to define the Gaudin model, in the case of periodic boundary 
conditions, by introducing the corresponding Lax matrix, 



L(A) 



N 

k=l 



(V.5) 



The Lax matrix obeys the Sklyanin linear bracket 



L(A),L( M ) 



ri 2 (A-f/), L(A) + L(p) 



(V.6) 
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A direct consequence of the Sklyanin linear bracket and the periodic boundary condi- 
tions is the fact that the operator 

t(A) = tr L(A) 2 (V.7) 

commutes for different values of the spectral parameter 

t(A)t(jO = t(^)t(A), (V.8) 

and is therefore the generating function of integrals of motion. The residues of t(A) at 
poles A = are the corresponding Gaudin Hamiltonians l6l lT6ll2T1| . 

A method to study open Gaudin model based on the following Lax matrix 

£(A) = L(A) - K(A) L(-A) K(-A), (V.9) 

and the corresponding r-matrix 

r K i 2 (A,f/) = r 12 (A- F )-K(A)r 12 (A + ^)K(-A), (V.10) 

was developed recently Il28ll45ll46ll . The relevant bracket in this case is given by 
f(A),£(p) 

The operator 

t(A) =tr£ 2 (A) (V.12) 

is considered to be the generating function of integrals of motion [28 45j[46]|. However, 
the Hamiltonians obtained as the residues of t(A) at the poles A = do not coincide 
with the Hamiltonians obtained in the previous section. 



r K 12 (A /f /),£(A) 



^2i(f/,A),£(» 



(V.ll) 



VI Conclusions 

We have reviewed the Gaudin model with integrable boundaries starting from the 
non-symmetric R-matrix of the spin- 1 /^ XXZ chain and generic solutions of the reflec- 
tion equation and the dual reflection equation. We have shown how the semi-classical 
expansion of the transfer matrix of the open inhomogeneous XXZ chain, calculated 
for special values of the spectral parameter, yields the Gaudin Hamiltonians with the 
boundary terms. Finally, an alternative approach to study open Gaudin model based 
on the classical Yang-Baxter equation and classical reflection equation is discussed. 
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